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Abstract. Many of the planets discovered via the radial velocity technique 
O \ are hot Jupiters in 3-5 day orbits with ~ 10% chance of transiting their par- 

^ • ent star. However, radial velocity surveys for extra-solar planets generally re- 

' quire substantial amounts of large telescope time in order to monitor a sufficient 

qq , number of stars due to the single-object capabilities of the spectrograph. A 

■ multi-object Doppler survey instrument has been developed which is based on 
the dispersed fixed-delay interferometer design. We present simulations of the 

, expected results from the Sloan Doppler survey based on calculated noise models 

■ and sensitivity for the instrument and the known distribution of exoplanetary 
^T) ' system parameters. We have developed code for automatically sifting and fit- 

. ting the planet candidates produced by the survey to allow for fast follow-up 

r***- ■ observations to be conducted. A transit ephemeris is automatically calculated 

C^l ' by the code for each candidate and updated when new data becomes available. 

The techniques presented here may be applied to a wide range of multi-object 

planet surveys. 

O ' 1- Introduction 

a: 

Of the transiting planets known thus far, three were discovered via the radial 
velocity t echnique and subsequently followed up ph otometrically. The se are HD 
209458b (Char bonneau et al. 1 120061 1 . HD 149026b dSato et al. II2005T ). and HD 



189733b (pouchy et al. 1120051 ), A multi-object Doppler survey instrument, the 
W.M. Keck Exoplane t Tracker, has been developed which builds upon the suc- 
<" pss ( Ge et al. 120061 ) of the prototype instrument, Exoplanet Tracker (ET), 



cess 

which is based on the dispersed fixed-delay interferometer design. This new 
instrument is being used to conduct an All Sky Extrasolar Planets Survey 
(ASEPS) with the Sloan 2.5m wide-field telescope, expected to dramatically 
increase the detection rate using the Doppler method. Since the initial survey is 
optimised towards the discovery of hot Jupiters orbiting stars in the magnitude 
range 8 < V < 12, the planet discoveries will be ideal candidates for the pos- 
sible detection of a transit sign ature in the corresponding lig htcurve. We have 



conducted survey simulations (jKane. Schneider, Ge 1 120061 ) to estimate how 



many transiting planets can be expected. 



Simulated Data 



To estimate how many planets we expect to detect in a given radial velocity ob- 
serving program, we performed a series of Monte-Carlo simulations which inject 
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Figure 1.: Cumulative planet-harbouring probability histogram (left) and the 
number of detectable planets as a function of RV precision (right). 



planets into a realistic sample of target stars based on the known distribution and 
characteristics of exoplanets. For this simulatio n, a stellar populat ion model was 
generated using the Besancon Galactic model ( Robin et al. 120031 ) for a magni- 
tude limited survey in the Kepler field. The Kepler field is an ideal location for a 
multi-object radial velocity survey since the results would not only compliment 
the transit survey to be undertaken by the Kepler mission, but also aid in target 
selection. In total, 23708 main sequence stars were used in the simulation from 
which 751 were concluded to harbo ur planets based on the planet-metallicity 
correlation ( Fischer Valenti 120051 ). The cumulative histogram of the planet- 
harbouring probability for the sample, shown in Figure 1 (left), demonstrates 
the substantial reduction in probability beyond ~ 5%. 

The distribution of planet parameters was used to calculate the expected 
distribution of radial velocity amplitudes. Figure 1 (right) shows the number 
of planets detectable from the simulated data at the 3a level as a function of 
the rms radial velocity precision of the experiment. In this example, doubling 
the precision of the instrument increases the planet yield by ~ 17%, whereas 
doubling the number of survey stars will increase the planet yield by 100%. 



3. Radial Velocity Fitting Code 

The FORTRAN code written for the purpose of sifting planet candidates from 
the dataset uses a weighted Lomb-Scargle (L-S) periodogram to detect a periodic 
signal in the data. If a significant periodic signal is detected, the data are 
subjected to an iterated grid-search fitting routine to determine the best fit 
planetary parameters for that target. Figure 2 (left) shows a typical periodogram 
where the dotted lines indicate various false-alarm probabilities. 

The number of false detections resulting from this technique depends upon 
the periodic false-alarm probability threshold one adopts as the selection criteria. 
Figure 2 (right) shows that by adopting a relatively conservative threshold, one 
can eliminate false detections with minimal loss to planet yield. 

Those planetary signatures failing to be recovered are generally at the pho- 
ton noise-limit of the instrument or have a period considerably more than the 
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Figure 2.: Periodogram for data containing a planetary signature (left) and the 
planets recovered as a function of the false-alarm probability threshold (right). 

survey duration. This sifting and fitting method has proven to be remarkably 
robust and is able to operate in an automated fashion. 

4. Number of Transiting Planets and Photometric Follow-up 

The total planet yield expected from the simulated data for the Kepler field sur- 
vey is ~ 275 planets for an instrument such as the ASEPS instrument assuming 
a 30 day observing window. To approximate the number of transiting planets, 
the geometric transit probability was calculated for each star/planet system as 
part of the Monte-Carlo simulation. Once all the detection limitations are taken 
into account, the total expected number is ~ 10 transiting planets. 

Considering the high number of expected hot Jupiters expected from the 
survey, the radial velocity fitting code automatically calculates a transit ephemeris 
for each candidate and is updated when new data becomes available. Since the 
transit duration is brief compared with the fitted period, the maximum win- 
dow for obtaining photometric transit observations after the radial velocity data 
has been obtained is also calculated. Through this process, we hope to provide 
a quick and efficient method for maximising the number of transit discoveries 
(lKanell2006h . 
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